function [cust_K] = custom_diffusion(xx,t)
    %return cust_K with the value of the diffusivity "K" for each node
    [nnod,ndm] = size(xx);
    cust_K = zeros(nnod,1);
    %% completar desde aqui
end